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We discuss the production of a class of heavy sterile neutrinos Vh in proto-neutron 
stars. The neutrinos, of mass around 50 MeV, have a negligible mixing with the 
active species but relatively large dimension-5 electromagnetic couplings. In partic¬ 
ular, a magnetic dipole moment ~ 10“® GeV“^ implies that they are thermally 
produced through e"'“e“ —>• i'h.i'h in the early phase of the core collapse, whereas a 
heavy-light transition moment Utr ~ 10“® GeV~^ allows their decay I'h with 

a lifetime around 10“^ s. This type of electromagnetic couplings has been recently 
proposed to explain the excess of electron-like events in baseline experiments. We 
show that the production and decay of these heavy neutrinos would transport energy 
from the central regions of the star to distances d ~ 400 km, providing a very efficient 
mechanism to enhance the supernova shock front and heat the material behind it. 
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I. INTRODUCTION 


Neutrinos define a sector of the Standard Model that still presents some important un¬ 
knowns. The current scheme of mass differences and mixings seems able to explain most 
of the existing data [l|, but the absolute value of their masses, their Dirac or Majorana 
nature or the presence of additional sterile modes j^, are yet to be determined. In 
particular, the production of sterile neutrinos Ug, through collisions with standard matter 
or flavor oscillations has important implications in particle physics jhl-Isl. The mixing with 
an active neutrino u may provide sterile modes with small couplings to the W and Z gauge 
bosons that translate into dimension-6 operators of type 


- ^es = - CAl^)f + h.c. 


( 1 ) 


In addition, the low-energy effective Lagrangian may also include dimension-5 operators 
from loops involving heavy particles. Although these operators are usually overlooked, they 
could mediate the dominant reactions of sterile neutrinos in a star under favorable ther¬ 
modynamical conditions. Here we will study this possibility in the context of supernova 
explosions. 

When a supernova goes off a proto-neutron star can be formed having typical initial 
radius (20-60) km and (1-1.5)Mq mass. It is believed that most of the gravitational binding 
energy (A’grav ~ 3 x 10^^ erg) is released in a ~20 second neutrino burst jo]. The neutrino 
spectrum from supernova SN1987A detected at SuperK and IMB indicated a weak decoupling 
from baryonic matter, conhrming that neutrino transparency sets in as their temperature 
falls below a few MeV [l^ in the dense core. At earlier phases of the collapse, however, 
computational simulations [ill, 1^ reveal internal peak temperatures exceeding 20 MeV in 
the central high density regions of the star. At such temperatures and densities the evolution 
of these astrophysical objects becomes sensitive to the fundamental properties of neutrinos 
and to the presence of hypothetical weakly-coupled particles. In this context, a lot of effort 
has been devoted to the cooling through neutrino emission in the nuclear medium [l^, to 
the matter opacity [l^ and revival of the stalled shock that arises in the standard paradigm 
of supernova core collapse 15[, or to the synthesis of heavy nuclei taking place in the hot 


bubble behind the shock M 


17 


In this work we will focus on the astrophysical consequences of the production of a heavy 
sterile neutrino Vh whose dominant interactions are not the weak ones in Eq. o but of 
electromagnetic kind. This type of particles have been proposed as a possible explanation 
18| for the excess of electron-like events in baseline experiments 19[. Let us briefly show 
how the required couplings could be generated. Consider a S'7/(2)2,-singlet Dirac neutrino, 
Vh, of mass m/j = 50 MeV. We will denote by N and the (2-component) neutrino and 
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FIG. 1. Diagram contributing to the magnetic dipole moment n of Vh. 
antineutrino spinors defining 



( 2 ) 


Let us also suppose that at TeV energies the gauge symmetry is SU{2)l x SU{2)ji x U{1)b-l 
and that Vh is accommodated within two SU{2)b, doublets together with a charged lepton, 


L = 




( 3 ) 


In order to avoid collider bounds 20(], the breaking of the left-right symmetry must be such 
that the charged lepton {E,E^) gets a mass rriE > 300 GeV while Uh remains light. Loop 
diagrams of heavy gauge bosons and fermions (see Fig. 1) will then generate the operator 


■^eff h" A , 


where is the electromagnetic field and /r is a magnetic dipole moment of order 211 


( 4 ) 
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( 5 ) 


In addition, the possible mixing of the sterile and the active neutrinos will be parametrized by 
an angle 6, so that the mass eigenstates read N' = cos 6 N+sin 6 v and v' = — sin 0 iV-l-cos 6 v. 
This mixing will generate electromagnetic transitions through the same type of diagrams (we 
drop the prime to indicate mass eigenstates): 


^ /itr Vh (1 - 75 ) u d^A’" + h.c., ( 6 ) 

with fitr ~ sin0/i being the transition dipole moment. This operator may imply that the 
dominant decay mode of the heavy neutrino is —)■ z/ 7 . Notice also that the presence of 
additional heavy singlets {vh' with mh' ^ mE) mixed both with u and Uh will give addi¬ 
tional contributions to /itr- Therefore, at this point we will treat /i and fitr as independent 
parameters. 








This type of sterile neutrinos could change substantially the evolution of a proto-neutron 
star. We will show that sterile pairs can be produced abundantly during the ~20 second 
neutrino burst, escape the star core more easily than standard neutrinos, and hnally decay 
within a few hundred km from the core. The very energetic photons from the decay could 
deposit energy, helping revive the stalled accretion shock formed during the collapse and 
change the thermal environment in the vicinity of the star. Our scenario could be considered 
a different realization of the eosphoric neutrino hypothesis proposed in 22 . 


II. DECAY RATE, PRODUCTION AND SCATTERING CROSS SECTIONS 


Let us first describe the dominant decay and production channels for the heavy neutrino 
Uh in vacuum. Later we will discuss how the hot and dense medium in a proto-neutron star 
(including populations of neutrons, protons, electrons and muons) affects these processes. 

To be dehnite in our calculation we will take as reference values rrih = 50 MeV, /i = 
10“® GeV“^ = 3.3 X and ptr ~ 10“® GeV“L We use h = c = 1. For these values of 

the mass and the transition moment, the heavy neutrino will decay into with a lifetime 

IOtt (50 MeV)3 (10“® GeV^)^ 

-o- ^ - 
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Fn ml 


mi 




X 0.0026 s . 


(7) 


We will also assume that the mass mixing, i.e., the active component in is smaller 
than sin^ 0 < 10“^ and only along the muon and/or the tan flavors. In that case, the 
radiative decay will dominate over the weak processes Uh —^ which appear 

with a branching ratio 


BR(i^/i —)■ ue'^e ) 
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( 8 ) 


10-3 (50 MeV)2 

This type of sterile neutrino avoids cosmological bounds since it decays before primordial 
nucleosynthesis. At colliders it is hardly detectable: even if it were produced in 1% of kaon 
or muon decays, Uh is too long lived to decay inside the detectors and too light to change 
signihcantly the kinematics of the decay 23|. Actually, more elaborate setups with two 
sterile modes have been proposed to explain the excess of electron-like events at MiniBooNE 
in terms of the photon that results from its decay 18 


The dominant production channels of will also be electromagnetic. In particular, 
electron-positron annihilation into Uh pairs, e-'-e” VhVh, will be mediated by a photon 
through the magnetic dipole moment coupling in Eq. (j4]). The differential cross section is 
given by 


da 

dt 
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FIG. 2. Left: Total cross section a{e~^e —)• Vh^h) for rrih = 50 MeV and /r = 10 ® GeV Right: 
a{uX —)• VhX) for fitr = 10“^ GeV“^ (solid: X = p, dashes: X = e). 

where a is the fine structure constant, mg is the electron mass, s and t are the usual Man¬ 
delstam variables is the center-of-mass energy). In Fig. 2 (left panel) we plot the total 
cross section for this process. Muon pair annihilation will give an analogous but subleading 
contribution, since muons are less abundant than electrons in the star core. 

The active to sterile transition mediated by a photon can be catalyzed by the presence 
of charged particles X = p,e: uX —)■ UhX (see right panel in Fig. 2). This contribution, 
however, can be neglected here due to the smaller value of the transition coupling that we 
have assurned, ptr ~ The weak channels which dominate the production of active 

neutrinos 


24| give also a subleading contribution due to the small mixing sin^ 0 < 10 ^ of 

25| are irrelevant for heavy neutrino 


our sterile, whereas other processes like plasmon decay 
masses around 50 MeV ( i.e., much larger than the electron mass). 

In addition to its production and decay, the collisions of with charged matter will 
be essential in order to understand its propagation in the dense medium and estimate how 
efficiently these neutrinos escape the proto-neutron star. We need to distinguish between 
elastic scatterings 

UhX UhX ( 10 ) 

and absorption reactions of type 

zz/jX —)■ z/X . (11) 

The differential cross section for the first process reads 
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FIG. 3. Left: cr{i'hX —>• for nih = 50 MeV, ^ = 10 ® GeV ^ and a scattering angle 0 > 30° 

in the center-of-mass frame. Right: a{vhX —)• vX) for fitr = 10“® GeV“^ (solid: X = p, dashes: 
X = e) 


This is a long distance (photon-mediated) process with a divergent total cross section; if we 
restrict to collisions substantially changing the direction of the incident { e.g., & scattering 
angle 6 > 30° in the center-of-mass frame) we obtain the cross section depicted in Fig. 3 
(left panel). For the inelastic process that transforms the heavy neutrino into an active one 
we obtain 
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dt 
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(13) 


and its total cross section is also depicted in Fig. 3 (right panel). These cross sections are, 
in both cases, much smaller than the ones for active neutrinos off nucleons mediated by 
weak bosons. For example, a neutrino in a gas at T ~ 20 MeV has an average energy 
{Ey) ~ itT ^ 60 MeV, and its lowest order elastic cross section with a neutron is a(vj?7, 
ViTi) ^ GjpEl{3C\ + Cy)/7T ~ 3 X cm^, with i = e, p,T (the cross section with a proton 
is approximately hundred times smaller). As for the absorption of a Vg through a charged 
current interaction, we have a{V(,n ep) 2 x 10“^® cm^. 


III. PRODUCTION IN A PROTO-NEUTRON STAR 


Let us now calculate the production rate of heavy neutrinos at the astrophysical site. 
Nucleon and lepton densities in the medium are constrained by electric charge neutrality 










and baryon and lepton number conservation. Typical baryonic densities in the core of a 
proto-neutron star are well above nuclear saturation density ub — (2-3)no with uq = 0.17 
fm“^ j^, whereas the lepton and electron fraction evolve dynamically from 0.31, 


Ye 0.27 at f = 0.1 s to Yl ~ 0.18, Ye ~ 0.17 at t = 10 s 26|. The extreme conditions in 


the core are such that quantum effects will be important. The population of baryons and 
leptons is described by the Fermi-Dirac distribution 


m) = 


^{E fJ.i)/T _|_ 2 ’ 


(14) 


where Hi {i = n,p, e^) denotes the chemical potential of the considered species. Both T and 
Pi evolve within the star, in particular, the chemical potentials take care of the conservation 
of charges and quantum numbers in a self consistent way 

For a given value of the temperature and the electron chemical potential {pe+ = “he-); fhe 
total energy emissivity (energy produced per unit volume and unit time) of heavy neutrino 

a 'rs through the dominant process —)■ Uhi^h-, QE{e'^e~ l^h^h) = is given by 


29 


Qe = 
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d^Pl Sp2 d^P3 d^P4 .4 


6 {pi + P2 - P3 - P4) (El + E 2 ) \M\ /(/l,/2,/3,/4) (15) 


(27r)8 J 2Ei 2 E 2 2 E 3 2 E 4 

where the factor /(/i, / 2 , /s, fi) = /i/ 2 (l — / 3)(1 ~ / 4 ) includes the Pauli blocking factor in 
the generic reaction 12 —)■ 34 and pi = {Ei,pi) are the 4-momenta. We will consider that the 
reaction is not affected by the quenching of outgoing sterile states, i.e., (1 —/s) ~ 1 ~ (1 —/ 4 ), 
since heavy neutrinos do not achieve chemical equilibrium and their number density inside 
the star is always small. The squared matrix element for the interaction dehned in Eq. (j4]) 
is given by 


lM{e^ 




-t 2 ml + ml 


— 2 {m\ + ml)t -I- {m\ — mlY 


(16) 


with e the electron charge, s = (pi + P 2 Y and t = {pi — Ps)^, see Fig. 4. 

To perform the phase space integral in Eq. flT^ we note that there are only four non¬ 
trivial independent variables: the initial energies Ei and E 2 , the angle 62 , as dehned in 
Fig. 4, and E 3 . Any other kinematical variables can be derived from these four or can be 
trivially integrated. It is convenient to dehne the 4-vector k = pi + p 2 (notice that k'^ = s 
and = (El -|- E 2 Y — s) and the angle 6*3 of k with p 3 \ 


2 (El - 1 - E 2 ) E 3 — s 

2\k\ ^E'i - mi 


(17) 


After integrating the 4-dimensional Dirac delta that enforces energy and momentum conser- 
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FIG. 4. Kinematical variables in the reaction 12 —)• 34 used in this work. 


vation, we obtain 
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^E+ 


P 1 P 2 


= / dcos02 / dEs^{E,+E2)\M\^feAEl)fe-{E2), 

64vr^ |fc| 

( 18 ) 

where the minimum values of E 2 and cos 6 * 2 , and c™™, respectively, result from the 

kinematical restriction s > 4m^ 


Er{E,) 
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The Mandelstam variables in terms of these four quantities read 


(19) 


( 20 ) 


s = 2 (ml + E 1 E 2 + P 1 P 2 cos 62 ) 
t = - 2 E 1 E 3 + 2 pip 3 cos(6*3 - a) +ml + ml 


( 21 ) 

( 22 ) 


where a (0 < a < tt) is the angle between k and pi, 

P 2 sin 6*2 


a = arctan 


Pi - P 2 cos 62 


(23) 


IV. TRANSPORT OF ENERGY OUT OF THE STAR CORE 


The possible impact of the heavy neutrino on the evolution of the proto-neutron star 
will depend on its ability to take a significant amount of energy out of the core. If the sterile 















neutrinos are abundant inside the star core but unable to reach the surface before decaying 
into a photon plus an active neutrino, then they become just a state mediating interactions 
of electrons with neutrinos and photons. We will show that this is not the case and that 
they could play an interesting role in supernova explosions. 

Let us take a temperature Tq = 25 MeV and an electron chemical potential = 100 
MeV, which are typical values at the inner central region of a proto-neutron star (see |lll.ll2[|). 
Although these quantities are time and density dependent, the chosen values can be used to 
estimate the possibilities of our scenario. 


Varying the mass ruh, the magnetic dipole moment /i, the electron chemical potential pe 
and the temperature T and performing a ht of Qe in Eq. (ll 8 l) we obtain 
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erg 
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(24) 


For the reference values of all the parameters, the expression above yields Qe ~ 2 x 
1035 gj-g^g cm^, with an average Uh energy of (Eh) 103 MeV. This is a very large pro¬ 
duction rate, ~ 10 ^- 10 ^ times larger than the one obtained in 22 | using heavier sterile 
neutrinos mixed with the active ones. Our neutrinos, however, will not leave the star core 
unscattered. 


We can also compare this production rate with the one of standard neutrinos in early 
cooling of proto-neutron stars. For example, in the central core the direct URCA process 
n —)> peue provides ^ 2.4 x lO^^i? erg/s cm^ at T = 25 MeV 28, M|, being R a 

factor of order unity [30|]. This is hve decades over the Uh production rate that we have 
found. The direct URCA process requires a high proton fraction, Yp > 11%, only accessible 
to large mass objects, however the less demanding modihed URCA cooling also gives a much 
faster rate than for steriles, ^ 1.5 x lO^^i? erg/s cm^. These active neutrinos will 

be, to a large extent, trapped (before transparency sets in) inside the star core, whereas ours 
have weaker interactions with the protons and electrons in the medium. 


It is then apparent that we need to consider propagation effects of Vh- Although a precise 
calculation would require a complete multidimensional simulation that is beyond the scope 
of this work, we will discuss the qualitative picture and show that the model has enough 
parameter space to realize it. 


The hrst important effect is the diffusion from the center to the surface of the star core, 
with a radius r ^ 20 km. Active muon and tan neutrinos scatter there mainly off neutrons 
with a mean free path ~ 1 m [3, . Analogously, will scatter elastically with protons 

with a cross section~40 times smaller, which suggests an interaction length Af inside the 
star longer by the same factor. This implies a larger diffusion coefficient D ^ ^hC/3 and a 
more efficient transport from the core to the outer parts of the star. The typical diffusion 












time for that process will be td ~ r‘^/{2D) 5 x 10“^ s. The second crucial effect in the 

propagation of the heavy neutrinos is their absorption: the star will capture a fraction of 
them through the inelastic collisions —)■ up. For /itr = 10“® GeV“^ the absorption length 

is approximately (j{uhX —)■ UhX)/a{uhX uX) 3000 times larger than Af, he., ~ 120 
km. A hnal effect, analogous to absorptions, is their decay. Since the values of ptv and rrih 
that we have assumed imply a sterile neutrino lifetime Th ~ 0.003 s, 16 times smaller than 
their diffusion time, a large fraction of the heavy neutrinos produced in the core will decay 
into yz/ before they have diffused to the outer layers. For a time window of ~ 20 s energy 
can be transported to distances d ~ \/2Dtd — 400 km. 

We estimate that, even if absorptions and decays reduced the number of neutrinos leaving 
the star from our estimate to 1% of the ones produced (some of them closer to the surface), 
Vh may still carry a total of lO^^-lO®^ erg and deposit this energy outside the star during 
the 20 second neutrino burst. Notice that a reduction in the magnetic moment p would 
also reduce the production rate of heavy neutrinos, but it would increase the mean free path 
between elastic scatterings and then the fraction of neutrinos that reach the surface. In 
addition to the transition moment ptr, the mass rrih is another parameter that could impact 
the production rate or the decay length of the heavy neutrino. 

An interesting variation of the model that depends only on two parameters {rrih and 
Ptr) would be obtained by suppressing the magnetic moment p and slightly increasing the 
transition one, e.g. /itr ~ 5 x 10“® GeV“^. In that case the dominant production channel 
would be 


uX —^ r'hX j uX —^ i^hX , 


(25) 


where X is any charged particle and the active u may be any linear combination of and 
Ut-. In the star core the heavy neutrinos would be partially absorbed through the inverse 
reaction (we estimate an absorption length Aj^ ~ 4 km), and the ones reaching the surface 
would decay with cTh ~ 30 km. Since the production would not be so abundant as in the 
case discussed in the previous section, this possibility provides a safer scenario still able to 
transport energy to the region near the star surface. 

Although the optimal value of these parameters would require a full Monte Garlo simula¬ 
tion, the scenario seems flexible enough to introduce acceptable changes in the dynamics of 
supernova explosions, with the decay into photon plus active neutrino playing an important 
role in the enhancement of the supernova shock. 



V. CONCLUSIONS 


The collapse of a very massive star defines an astrophysical object with extreme conditions 
where neutrinos determine the thermodynamics. These proto-neutron stars are suitable 
laboratories to probe the properties of any weakly coupled particles of mass ^ 100 MeV. 
Here we have proposed a sterile neutrino much heavier than the standard ones {rrih ~ 50 
MeV) and with sizable electromagnetic couplings: a magnetic dipole moment p and dipole 
sterile-active transition fitr that mediates its decay Uh —t zzy with a lifetime ~ 10“^ s. This 
simple 3-parameter model seems to have interesting implications in supernova explosions. 
The heavy neutrino is produced in the core at a high rate through Vh^hi h may 

escape the star more efficiently than active neutrinos and decays depositing a large amount 
of energy in the outer layers of the star. 

We believe that the type of heavy sterile neutrino proposed here could be an essential 
ingredient to help the progression of the internal shock, which is responsible for the observed 
supernova events. Full computational simulations could shed more light into the complex 
energy transport that results from competing processes of scattering, interaction and decay. 
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